Quantum phases of a multimode bosonic field coupled to flat electronic bands 
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We investigate the quantum phases of systems in which a multimode bosonic field is coupled to 
the transitions between two flat electronic bands. In the literature, such systems are usually modeled 
using a single or multimode Dicke model, leading to the prediction of superradiant quantum phase 
transitions for large enough couplings. We show that the physics of these systems is remarkably 
richer than previously expected, with the system continuously interpolating between a Dicke model 
exhibiting a superradiant quantum phase transition and a quantum Rabi model with no phase 
transition. 



In cavity quantum electrodynamics experiments the 
photon confinement increases the light-matter coupling, 
eventually allowing for a single photon to be absorbed 
and re-emitted multiple times before escaping out of the 
cavity. It is the hallmark of the strong coupling regime, 
realized for the first time with Rydberg atoms in super- 
conducting cavities Still, while the strong coupling 
regime is observable in these systems thanks to the im- 
pressive quality factors of superconducting cavities, the 
strength of the dipolar coupling between a single Ryd- 
berg atom and the electromagnetic field, quantified by 
its quantum Rabi frequency, is bounded by funda- 
mental reasons to be a small fraction of the frequency of 
the transition, CJ12 [2j. 

This bound can be lifted if photons couple to collec- 
tive superradiant electronic excitations [3]: TV identical, 
coherently excited dipoles, behave as a single super dipole 
y/~N times larger. Such a phenomenon opens to the possi- 
bility to experimentally achieve the ultrastrong coupling 
regime, in which the vacuum Rabi frequency becomes of 
the same order of the bare frequency of the electronic 
excitations, and the normalized coupling, is thus of 
order one [J-l9| . As in this regime the interaction energy 
is one of the dominant energy scales of the problem, it 
becomes possible to observe fascinating new physical phe- 
nomena, ranging from quantum vacuum radiation JJJ- 1_2| 



to superradiant quantum phase transitions (QPT)[13l- 

These systems are usually studied using the Dicke 
model [U, describing TV two level atoms identically cou- 
pled to the electromagnetic field. This model is based on 
the assumption that the interaction with the photonic 
field can only induce transitions between levels in the 
same atom. No transitions between different atoms are 
allowed. While this assumption is well justified for dilute 
atomic clouds, this may not be the case in solid state 
systems @, [?J or in synthetic many-body systems [T9I — 
[22[, where the two level systems can form, instead, two 
electronic bands. In this case, each state of the lower 
band can be coupled, through different photonic modes, 
to different final states in the upper band. 

The Dicke model can only give an approximate de- 
scription of these band models, as it can be grasped by 



a simple dimensional argument. The Hilbert space of 
TV electrons in a Dicke model has dimension g^dm = 2^ 
(TV two level systems). For the same TV, the electronic 
Hilbert space in a band model, in which electrons can 
jump from any state of the lower band to any state of 

' 2 TV ' 



the upper one, is instead c^fbm = 



TV 



(TV electrons 



spread over 2TV possible states), that tends to 



tN 



for 



large TV. The physical Hilbert space is therefore much 
larger than the one of the Dicke model and, as we will 
see, this translates in a much richer physics. 

In this Letter, in order to investigate this new and com- 
plex physics, we will study a model describing a multi- 
mode bosonic field coupled to two flat electronic bands, 
separated by a bandgap lj 12 • This particular choice is 
motivated by the fact that flat band models, that can 
be engineered using a number of different systems, from 
two dimensional electron gases [7|, |8| and graphene under 



magnetic field [23|, [24| , to atoms in optical lattices [2ol — 
[22| and topological insulators 25l-[28|. are characterized 
by density of states peaked around the badgap, leading 
to the possibility to observe extremely large light-matter 
couplings [111, [23|. 

We will prove that, modifying the spectrum of the 
bosonic field, it is possible to continuously tune the 
physics from the usually expected Dicke model to an ul- 
trastrongly coupled quantum Rabi model [29l-[32|. Hence, 
the engineering of the bosonic spectrum can allow to sim- 
ulate different classes of superradiant models, with a dra- 
matic impact on the phase diagram. 

We consider two flat bands of states with TV modes 
each, with the Fermi level in the bandgap. The states in 
each band are indexed by a general n-dimensional mo- 
mentum, k, with k £ [— fco,fco], and periodic boundary 
conditions. Each state in the lower band, c^k, is cou- 
pled to each state of the upper band, C2,k', by the ap- 
propriate momentum-conserving bosonic mode, a q , with 
q = k 7 — k, through a coupling constant, f2, independent 
from q (see Fig. [T] (a) for a schematic representation of 
the model). Such a model is of course an idealization, and 
the physics of a real implementation could be modified by 
the presence of cut-offs and inhomogeneities, due both to 



2 



the Fermi surface and to the wavevector dependency of 
the coupling constant. Still, this idealized model allows 
us to show, in a clear and paradigmatic way, how the 
presence of extra electronic transitions, coupled to dif- 
ferent modes of the bosonic field, can profoundly modify 
the physics with respect to the well-known Dicke model. 




FIG. 1: (a) Sketch of the considered model with flat bands 
in momentum space: electron can occupy two parallel bands, 
and each state in the first band couples to every state of the 
second one through a different bosonic mode, with a coupling 
constant -^=. (b) Real space picture: each electron in the 
first band couples with only one state in the second subband, 
through a local bosonic field, with a coupling strenght Q. The 
different local bosonic modes at sites 1*1 and r2 are further 
coupled between them by the coupling constant Cj r2 -r 1 - 



The Hamiltonian of the system, in momentum space, 
thus reads 

H FBM = UqQ ^ + ^2C\ M C 2M (1) 

q k 

+ TT^f X^ C 2,k+q C l,k + ct l,k C 2,k-q)(a t -q + Gq), 
k,q 

where uo q is the frequency of the mode a q , that we will 
take to be isotropic, and we have set, here and in the 
rest of this paper, h = 1. As we have anticipated in the 
introduction, the behavior of the Hamiltonian in Eq. (pQ) 
is controlled by the dispersion of the bosonic field, uo q . It 
is thus practical to introduce three quantities describing 
its average, oj = ^] q a; g , its half width, A = uo — cj m , 
and the energy of the lower lying mode, uo m = min q uo q . 

Defining the Fourier transform h Yl of a quantity /i q , 
as hq = ^= J2 r ^re 2qr (r is thus a point on an n-cubic 
lattice of side |^), we can rewrite the Hamiltonian in Eq. 



(PP) in real space as 

H FBM = J2u4a r +LJ 12 ^ r C2,r (2) 
r 

+ 0(4 r^l,r + c\ r C 2 , r )(4 + CL r ) + V] ^^a\a Y > , 
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where we have exploited the fact that uo = -^=. 

Restricting the Hilbert space to states in which only 
one electron is present at each site r, we can map the 
matter part at each site to a two level system. This is 
equivalent to neglect pairs of sites, one empty and one 
doubly occupied, completely decoupled from the electro- 
magnetic field. This simplification can be justified, a 
posteriori, by the fact that the energy of such a pair, 
GJ12, is much higher than the energies we will find for two 
interacting, singly occupied sites. 

Hamiltonian in Eq. (j2j) is reminiscent of a Rabi- 
Hubbard model [111 , describing a lattice of two- level sys- 
tems, each coupled to a local bosonic mode of frequency 
lj, with a superradiantly enhanced coupling strength ft. 
The coupling of two sites at distance r being given by the 
r-component of the Fourier transform of the boson dis- 
persion uo r (see Fig. [T] (b) for a schematic representation 
of the model in real space). 

For a flat bosonic dispersion, A = 0, uo q = uj m = a), and 
the only nonvanishing Fourier component corresponds to 
r = 0. Different sites thus completely decouple and in 
this limit the Hamiltonian in Eq. (j2j) exactly maps on the 
sum of N identical and decoupled quantum Rabi models 

H QRM = u m rfa + ^y(j z + ^(a f +5), (3) 

where <r J , j G [x,y, z] are Pauli matrices. Notice that in 
Eq. ([3]), and in all the rest of this Letter, we discard the 
constant terms from the Hamiltonian and neglect the site 
index for fully local and site-independent Hamiltonians. 

This exact mapping has two, far-reaching conse- 
quences. On one side, even if Eq. (j3j) describes the 
physics of a single two level system, the light-matter cou- 
pling is enhanced by the superradiant factor y/N . Using 
flat bosonic dispersions, like the ones that can be engi- 
neered with photonic crystals, a flat band model can thus 
be used to simulate an ultrastrongly coupled quantum 
Rabi model, opening to the possibility to observe new 
physics, as quantum vacuum radiation [34] or the Bloch- 
Siegert shift [35]. On the other side, while the Dicke 
model presents a phase transiton, when the coupling ft 
exceeds a certain critical value [H, Qjjl, EH , the quantum 
Rabi model does not. We thus proved that the Hamilto- 
nian in Eq. (pp) is in not always in the same universality 
class of the Dicke model. 

In the non-flat dispersion case however, the phase tran- 
sition can occur. In order to prove this explicitly, we will 
use a mean-field approach to study the full Hamiltonian 



3 



in Eq. (j2j). We shift the bosonic fields as 

aJ = v>r + <£, (aj) = ^r, (<£) = o, (4) 

where the expectation value is taken over the ground 
state. In order to decouple the different sites, we disre- 
gard the coupling between fluctuations at different sites 
and obtain, on each site r, the following Hamiltonian 



MF 



^12 



(5) 



We expect the previous mean field Hamiltonian to 
present a delocalized superradiant phase transition [33j| 
for large enough values of the inter-site coupling. In order 
to allow for the field to condense in any of the original q 
modes, we chose the following inhomogeneous ansatz for 
the mean field parameter 



= ^>/2cos(qr) for q ^ 
^ = if) for q = 0, 



(6) 



with tp real and the g-dependent normalization is needed 
to have a continuous energy for the bosonic field. 

Inserting the ansatz in Eq. (|6| into Eq. (j5| we obtain 



H, 



MF 



^12 

2 



(7) 
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FIG. 2: Phase diagram of the considered flat band model 
at the resonance cj m = CJ12, as a function of the normalized 
coupling, , and of the normalized half width of the bosonic 

dispersion, We see how the normalized coupling at the 

critical point is 0.5 in the limit of wide bands, as predicted 
by the Dicke model and diverges for flat bands, as the system 
becomes well described by a quantum Rabi model. 



In order to calculate the phase diagram of the system 
we have to find the values of ip and q that minimize the 
ground state energy of the total mean field Hamiltonian 
^fbm = Xlr '• As a first step we notice that the 
Hamiltonian in Eq. ([7]) is invariant over a contemporary 
sign exchange of a* and <5j. The expectation value 
of an arbitrary observable Tfr, over site r, can thus be 
written, close enough to the transition to be able to keep 
only the lowest order term in the order parameter, in the 
form 



q2 



(8) 



where A q and B q are independent of ip^ (but they can 
still depend upon q through uj q ). An homogeneous ob- 
servable, 77, is equal to the average, over all sites r, of 
local observables r] r . From Eqs. (|6j) and (|5J), we thus 
have that 



7] = A q + Bqi/j 



(9) 



For the sake of determining the phase boundary, we can 
thus consider the homogeneous case, dropping the site 
index. Moreover, if we calculate the ground state expec- 
tation value, we obtain from Eq. (J2J) 



(H Mb ) 



,2 , / w 12 



na x {5^ + 5)). 



(10) 



The first term in Eq. ([TO]) is a strictly increasing function 
of uo qi while the others do not depend upon q. The sys- 
tem will thus always condense in the lowest lying bosonic 
mode and, now on, we can take uo q = uo m . We can thus 
limit ourselves to study a single, site independent local 
Hamiltonian 



H 



MF 



^12 



(11) 



An analytic limit for the phase boundary can be ob- 
tained for broad dispersions (such that the average en- 
ergy mismatch is much larger than the coupling, A ^> f2). 
In this limit, we see from the fourth term in Eq. ([TT]) 
that the energy cost of fluctuations increases, eventually 
completely freezing them. It is thus possible to neglect 
quantum fluctuations and consider only the coherent part 
of the field 



lim H 



W12 



-^cj z +cj m ^ +2ta/<. (12) 



The previous Hamiltonian can be analytically diagonal- 
ized, yielding an energy, close to the QPT boundary 
E = (uj m — ^-)V ;2 7 leading to a second order QPT for 
the critical value 



v /o; m cJi2 



(13) 



that is the critical value we expect from the Dicke model 
13l , [HI- This result justifies, a posteriori, the choice 
to model flat band systems with steep enough bosonic 
dispersions as Dicke models. In this case, all the bosonic 
modes except the almost resonant one eventually become 
completely nonresonant and can thus be neglected. 

For intermediate values of A, we diagonalized H MF in 
a truncated Fock space for different values of finding 
then the value that minimizes the ground state energy 
36|, [37j . In Fig. [2]we plot the phase diagram as a function 
of the normalized coupling, and of the normalized 
half width of the bosonic dispersion, ^ , for the resonant 

case uo m = uj\2. In the two limiting cases of going to- 
ward zero and infinity, we find results consistent with 
our previous analytical estimations: the critical point di- 
verges toward infinite values of the normalized coupling 
when the dispersion is almost flat (consistently with the 
fact that the quantum Rabi model does not present a 
critical point) or it converges to 0.5, accordingly to Eq. 
(fT3|) . for very broad dispersions. 

In Fig. [3]we show the average occupation of the excited 
band, X^(4,k c 2,k) = (cr z ) + |, for the same param- 
eters. From the bottom panel of Fig. [3] it is evident 
the continuous transition between the Dicke-like quan- 
tum phase transition for broad dispersions (solid line) 
and the Rabi-like crossover (dotted line). 

In conclusion, we investigated the physics of a model 
describing a multimode bosonic field coupled to two flat 
electronic bands. We proved that, if the bosonic mode 
is also dispersionless, the presented model maps exactly 
into an ultrastrongly coupled quantum Rabi model, al- 
lowing to simulate such a model with controllable cou- 
pling strengths. Using a mean field theory, we showed 
that, engineering the dispersion of the bosonic field, the 
physics of such a model can be continuously tuned from 
a quantum Rabi model (with no critical point) to a Dicke 
model exhibiting a superradiant QPT. 

We thank D. Hagenmiiller, J. Keeling and P. Nataf for 
useful discussions. 
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FIG. 3: Average occupation of the states in the excited band, 
for ujm = CJ12. In the top panel we show a colorplot as a func- 
tion of the normalized coupling, and of the normalized 
half width of the bosonic dispersion, In the lower panel 

we show instead four cuts of the same image, along the four 
vertical lines visible in the colorplot. We see how the system 
interpolates between the continuous behavior of the quantum 
Rabi model and a second order quantum phase transition for 
the Dicke model. 
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